Magnetic Vortices in High Temperature Superconductors 
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It is suggested that modes, observed in recent neutron scattering experiments by Lake et al, 
on La2-a;SriCu04 in strong magnetic fields (« 7 T), are due to the existence of antiferromagnetic 
moments associated with the cores of vortices generated by the field. These moments form one- 
dimensional chains along the c-axis (the vortex axis), which at finite temperatures are disordered. 
At temperatures higher than 10 K the correlation length gets shorter than the lattice parameter, 
resulting in no scattering from coherent spin-waves above that temperature. The bandwidth of the 
spin-waves is estimated to be « 4 meV in accordance with the observations. 
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The superconducting and the antiferromagnetic phases 
are close in the copper oxide supersuperconductors. They 
are close in the overall phasediagram, where at low tem- 
peratures the antiferromagnetic (AF) order is replaced by 
a superconducting (SC) order upon doping and in certain 
inhomogenous phases, such as the striped phases, the two 
order parameters will be literally close geometrically. In 
the original observations of stripes in cuprates by Tran- 
quada et al. it was found that the stripes are in fact 
insulating regions with antiferromagnetic order in a sea 
of superconducting material. An alternative way of sup- 
pressing the superconductivity is to apply an external 
magnetic field which will generate vortices, who's cores 
are non-superconducting. It has been suggested by sev- 
eral, that there is a possibility, that the cores are in fact 
small islands of insulating antiferromagnets. 

So far we have only indirect evidence, that such a pos- 
sibility is realized. E.g. Scanning Tunneling Microscope 
(STM) measurements Q indicates an insulating behavior 
of vortex cores in the Bismuth based cuprates. Recent 
inelastic neutron scattering experiments m on optimally 
doped La2_3;Sr3;Cu04 in an applied magnetic field has 
revealed modes with energies around 3.5 meV below 10 
K. In this paper it will be argued that such results are to 
be expected if in fact the vortex cores support antiferro- 
magnetically ordered spins. 

In the experiments by Lake et al. on optimally 
doped La2-2;Sr2;Cu04 below Tc it is found that with- 
out a field the spin excitations have a gap A = 6.7 meV, 
for all measured momentum transfers. A simple BCS 
model for c?-wave superconductivity 1^, leads one to ex- 
pect that the gap in the spin-excitations should vanish for 
momentum transfers close to the difference between two 
nodes in the one particle spectrum. This is not found ex- 
perimentally, from which one concludes, that in strongly 
correlated systems the spin excitations are not necessar- 
ily associated with the charge excitations. We shall have 
nothing further to say about this general problem. When 
a field is applied, however, the situation changes. Then 
new excitations are found below a temperature Tq = 10 
K. These excitations have energies in the gap associated 



with the field free superconductor. In fact the results can 
be summarized by saying that there is a single mode with 
energy 3.5 meV and width 3.8 meV, although this is not 
the physical picture emerging from the considerations in 
this paper. It is natural to associate the excitations with 
the vortices generated by the magnetic field, and Lake 
et al. also establishes that the intensity of the scattering 
scales with the strength of the external field, consistent 
with such an assignment. 

Using the generalized Ginzburg-Landau model, which 
includes both the AF order parameter and the SC order- 
parameter, and which has an approximate SO (5) symme- 
try, Arovas et al. Q has suggested the that vortex cores 
can be insulting antiferromagnets. The scenario is not 
only possible in an S0(5) model, and more recently it 
has also been proposed in models with spin charge sepa- 
ration. In this paper we shall make a minimal assump- 
tion, that in the vortex core the AF order parameter is 
non-zero in the core, and described by an envelope func- 
tion 4>{r), who's total weight J (t>{rYrdr / = TV, where 
N then can be interpreted as the total number of spins in 
the vortex core. The excitation of such an AF vortex core 
has been studied in the work by Hallundbaek et al. 
Neglecting for a moment the interaction between spins in 
different copper-oxide layers and spin anisotropics, the 
spins in the core has a zero-energy, Goldstone mode, 
which has excitation energy zero, corresponding to a rigid 
rotation of all the spins in the core. This mode is model 
independent. Since the group of spins has a finite size, 
the magnons will have a discrete spectrum, and the gap 
to the lowest excitation is approximately Ja/^, where J 
is the exchange interaction in the plane (~ I'iQmeV)^ a 
is the lattice parameter, and ^ is the coherence length, 
i.e. the size of the vortex core, hence the energy is much 
larger than the few meV we are considering in this pa- 
per, and we can take the spins to form a rigid group 
described by one single vector n. Different models have 
different other excitations of the core. E.g. will the S0(5) 
model have a resonance close to 40 meV, which describes 
a rotation of the AF order parameter into the SC plane. 
The details are not important — here we only need to 
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know that the energy is so high, that such modes can 
be neglected in the present context. There are two final 
interactions, which have low energies: the Zemann cou- 
pling to the magnetic field associated with the vortex, 
and the interplane exchange coupling J'. Knowing the 
London length one can easily estimate the Zemann en- 
ergy, and this is 1 /ieV, which is way to small to be of 
relevance This leaves us with J'. It is known to be 
J' = 3.7 • 10~^ J, which also seems small, but it is an in- 
teresting feature of spin waves in a system with different 
exchange constants in different spatial directions, that 
modes moving along the direction with J' (the c-direction 
in our case) have energies ujk = \/ J J' {I — cos fee), i.e. a 
typical energy V J J' « 2meV, which is the energy scale 
of the experiment. 

So, we are suggesting, that what is seen in the ex- 
periment is spin waves traveling up and down the vor- 
tex core column. From general principles we know that 
a 1-dimensional antiferromagnet does not order at any 
finite temperature, so we have to analyze a thermally 
disordered antiferromagnet. This is most easily done us- 
ing the field theory representation for an antiferromag- 
net, namely the non-linear sigma model. We start with 
the model in the 2-dimensional copper-oxide planes. The 
(imaginary time) action is given by 



1 / 1 / dm 



2\AaJ \ dr 



+ aJs^{V'mi)^ +Lsc\ , (1) 



with X = N/{4aJ) and p = J's'^cN/A. 

Both the "mass" and the "spring constant" scales with 
the size of the vortex core. A'', but at zero temperature 
this model have low energy spin wave modes with a dis- 
persion iujn = Sy^ JJ'a/c kc, independent of A'^ — and 
equal to the dispersion of spin waves in an 3-D model 
with anisotropic exchange constant. The size of the vor- 
tex core can be estimated from recent measurements of 
Hc2 by Ando et al. § They find Hc2 ~ 62T, which result 
in an estimate for A^ « 216. 

At finite temperatures the system is disordered. The 
basic physical picture put forward in this paper, is that 
only propagating modes with a wavelength shorter than 
the correlation length will show up as peaks in inelastic 
neutron scattering. Below a certain temperature (10 K 
in the experiment) the correlation length becomes suf- 
ficiently long, that modes with wave vectors in the ex- 
perimental window (« 50 % of the c-axis Brillouin zone) 
contribute to the scattering. Such modes has energies in 
the upper half of the spin wave spectrum. 

Accordingly we need to find the correlation length of 
the spins. We use the renormalization group equations 
discussed extensively by Chakravarty et al. Q . They are 
for the 1-D case 



§ = f cothff) 
dl 2tt \2tJ 



dt 



2n \2tJ 



(5) 
(6) 



where g is a dimensionless coupling constant: 



where s is the size of the spins, a is the in-plane lat- 
tice constant and i is a plane index. The last term 
describes the coupling to the SC orderparameter. The 
precise nature of that is not specified, except, that we 
assume, that in the core of the vortex, the AF order 
parameter TOi(r) is a constant vector n;, i.e. we make 
the ansatz mi(r) = ni0(r), where 4>{r) is 1 in the cen- 
ter of the core, and decays to zero outside the core. 
The actual expectation value of the spin-operator is then 
{Si{Rj)) = sni(t>{Rj){-iy , where the factor (-1)^ de- 
notes the AF staggering, +1 on one sublattice and -1 on 
the other sublattice. 

Let us now introduce the inter plane coupling: 
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and t is a dimensionless temperature: 



t 



NJ's"' 



(7) 



(8) 



The correlation length for the spin system, is deter- 
mined by the requirement that the renormalized length, 
^s(0 — ^sG~^ is equal to the lattice parameter c This 
is essentially obtained when the scaling parameter, / has 
a value such that the renormalized temperature is 27r^. 

The renormalization equations can be solved exactly, 
and the solution is 



which can be rewritten in terms of the in-plane AF order 
parameter (and making the continuum approximation): 



J's^cN 



dz 



4 J V^^. 
resulting finally in the effective action 

I r'^ . f . ( fdnV fdr 




(3) 
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where 



cothhf^e-'' 
27r \2to 



(9) 



The correlation length is finite at all temperatures, with 
a limiting value 



= = r.p27r/go 



UT = 0)=ce 



(10) 
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For the parameters relevant to our system, the scahng 
starts at values 

50 ~ 5 (11) 
io«1.7T, (12) 

where T is the actual temperature measured in Kelvin. 
These are approximate values derived from very difBcult 
experiments (iJc2 and J' are hard to get), so precise scale 
for e. g. to can easily vary with a factor of 2. Chosing the 
temperature scale factor to bo 3.4 instead of 1.7, which 
corresponds to the product NJ' being a factor 2 smaller 
than quoted above, we got the curves shown in Figure 1 
for the correlation length ^(T) as a function of temper- 
ature and for 3 values, 1, 5, 10 for ^o- We see that a 
strong growth of correlation length around 10 Kelvin is 
quite compatible with known values for the size of the 
vortex core and the strength of the interplane exchange 
coupling. 



approximately 3 mcV. This energy scale is the natural 
one for spin waves in a system of antiferromagnetically 
ordered spins residing in the vortex cores of the supercon- 
ductor. This being a 1-dimensional system no long range 
order exist at any finite temperature, so inelastic scat- 
tering creating reasonably longlived spin waves requires 
a correlation length larger than the interplane distance. 
We have shown that for reasonable values of the vortex 
core size and the interplane exchange coupling, this hap- 
pens for temperatures less than 10 Kelvin. It is beyond 
the scope of this paper to give a detailed calculation of 
the actual theoretical neutron scattering spectra, but cal- 
culations of this is underway. 

I thank Bella Lake and Kim Lefmann for an early peek 
at their data. 
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FIG. 1. The correlation length of spins in the vortex core 

along the c-axis as a function of temperature. The lower curve 
is for go — 10, the middle for go = 5 and the upper for qq = 1. 
The temperatur scale is given by to = 3.4T, where T is mea- 
sured in Kelvin. 

In conclusion, we have proposed an explanation of 
the recent experimental finding by Lake et al. of mag- 
netic scattering from optimally doped La2-xSrxCu04 in 
a strong magnetic field, below 10 Kelvin at an energy of 
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